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Abstract 

This paper describes perturbative framework, on the basis of the closed- 
time-path formalism, in terms of quasiparticle picture for studying quasi- 
uniform relativistic quantum field systems near equilibrium and nonequilib- 
rium quasistationary systems. Two calculational schemes are introduced, 
the one is formulated on the basis of the initial-particle distribution function 
and the one is formulated on the basis of the "physical" -particle distribution 
function. It is shown that both schemes are equivalent and lead to a gen- 
eralized kinetic or Boltzmann equation. Concrete procedure of computing a 
generic amplitude is presented. 
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1 Introduction 



Since mid-fifties, efforts have been made to incorporate quantum field theory with 
nonequilibrium statistical mechanics]!], |2j, [3|, ||. Necessity of this incorporation origi- 
nated from the field of solid-state physics. Since then, rapid progress of the studies of 
the early universe and the quark-gluon plasma, which is expected to be produced in 
heavy-ion collisions and to have existed in the early universe, have further activated this 
field of research (cf., e.g., Refs. 5) and 6)). 

The standard formalism of nonequilibrium statistical quantum-field theory is broadly 
classified into two frameworks, the one is the closed-time-path (CTP) formalism [III [2|. 
|3], [7|, and the one is nonequilibrium thermo field dynamics.^] In this paper, we employ 
the former. Throughout this paper, we are interested in quasiuniform systems near 
equilibrium or nonequilibrium quasistationary systems, which we simply refer to as out- 
of-equilibrium systems. A brief discussion is made in §5 on the relation with thermo 
field dynamics. 

The out-of-equilibrium systems are characterized by two different spacetime scales: 
microscopic or quantum-field-theoretical and macroscopic or statistical. The first scale, 
the microscopic-correlation scale, characterizes the range of radiative correction to re- 
actions taking place in the system while the second scale measures the relaxation of 
the system. For a weak-coupling theory, in which we are interested in this paper, the 
former scale is much smaller than the latter scale. A well-known intuitive picture (cf., 
e.g., Ref. 7) for dealing with such systems is to separate spacetime into many "cells" 
whose characteristic size, (// = 0,...,3), is in between the microscopic and macro- 
scopic scales. It is assumed that the correlation between different cells is negligible in 
the sense that microscopic or elementary reactions can be regarded as taking place in a 
single cell. On the other hand, in a single cell, relaxation phenomena are negligible. 

Above intuitive picture may be implemented as follows. Let A(x,y) be a generic 

propagator. For an out-of-equilibrium system, A(x,y), with x — y fixed, does not change 

appreciably in the region \X^ — Xq | < L M , where = (x^ + y fl )/2 is the mid-point and 

Xq is an arbitrary spacetime point. A self-energy part H(x,y) enjoys similar property. 

Thus, X^ may be used as a label for the spacetime cells and is called the macroscopic 

^The nonequilibrium thermo field dynamics is initiated by Arimitsu and Umezawa^j and has devel- 
oped into sophisticated form. j|, [n| 
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spacetime coordinates. On the other hand, relative spacetime coordinates x^ — are 
responsible for describing microscopic reactions taking place in a single spacetime cell. 
An inverse Fourier transformation with respect to the relative coordinates x — y yields 



A(X;P) = J d\x-y)e iP <^A(x,y) (1.1) 



(with P M = (p°,p)) together with a similar formula for E. Above observation shows 
that, for |P M | > l/L^, P M in (|1 . 1|) is regarded as the momentum of the quasiparticle 



participating in the microscopic reaction under consideration. 

Microscopic reactions discussed above cause change in the density matrix that charac- 
terizes the ensemble of the system, through which the number densities of quasiparticles 
change with macroscopic spacetime A M . Dealing with this is the subject of the "next 
stage," where (weak) A M -dependence of A(X; P) and T,(X; P) are explicitly taken into 
account. 

The purpose of this paper is to clarify the structure of perturbative out-of-equilib- 
rium quantum-field theory through introducing two mutually equivalent calculational 
schemes. The one (bare-iV scheme) is formulated in terms of the initial- or "bare"- 
number-density function (§2 - §3) and the one (physical- N scheme) is formulated in terms 
of the "renormalized" -number-density function (§4). We demonstrate the equivalence of 
both schemes through deriving a generalized Boltzmann equation. Concrete procedure 
of computing a generic amplitude is presented in terms of Feynman rules supplemented 
with the derivative expansion and the generalized Boltzmann equation (§4). The bare- 
N scheme is the scheme which is directly deduced from first principles. To the best 
knowledge of the present author, this scheme has never been comprehensively dealt 
with. The physical- N scheme is the scheme which is employed in the literature. Within 
the CTP formalism, however, no self-consistent deduction of the scheme is available. We 
deduce this scheme from first principles in an unambiguous manner. Comparison with 
related work is made in §5. 

Throughout this paper, for the sake of definiteness, we take up a self-interacting 
complex-scalar quantum field theory. Generalization to other theories is straightforward. 
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2 Closed-time-path formalism 



To make the paper self-contained, in this section, we give a brief review of the CTP 
formalism. 

2.1 Preliminaries 

The Lagrangian (density) reads 

£ = -0t(9 2 + m 2 )0-^(0t0)2. (2.2) 

The CTP formalism is constructed through introducing an oriented closed-time path C 
(= C\®C2) in a complex-time plane, which goes from — oo to +oo (Ci) and then returns 
back from +oo to — oo (C2). The time argument Xq of the fields, <f>(x), is on the time 
path C. The perturbative scheme for computing time-path-ordered Green functions is 
summarized as a set of Feynman rules (cf. Ref. 2)). 

Throughout this paper, we assume, for simplicity, that the density matrix p commutes 
with the charge operator Q: [p, Q] = 0. (Generalization to the case [p, Q] 7^ is 
straightforward.) Building blocks of Feynman rules are the propagator, the vertex and 
the so-called initial-correlations. The propagator A is defined by 

A(x,y) = -iTr [T c (<f>(x)^(y)) p] (= -i(T c (<J>(x)4>\y)))) (2.3) 

and the initial correlations are 

C 2 n = i(-) n Tr[:^(x 1 )... ( f ) (x n )^(y 1 )... ( f ) \y n ): p] (n > 2). (2.4) 

Here T c is the time-ordering operator along the time-path C and : ... : indicates to 
take normal ordering with respect to the creation and annihilation operators in vacuum 
theory. Throughout this paper, we do not deal with C 2n (n > 2). (In Appendix A, we 
discuss a condition, under which C 2n (n > 2) may be ignored when compared to the 
2n-point Green function.) 

A field 4>(xq,x) with xq G C\ [C2] is called a type-1 [type-2] field, which we write 
<t>i{x) [02 (x)]. The type-1 field is also called a physical field. A classical contour action 
may be written in the form, 

j dx Q j d 3 x£(0(x),0t( x )) = J + °°dx J d^x^i^x),^)) - C{<t> 2 {x),<t>\{x))\ 
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d 4 xC(x). (2.5) 

Following Refs. 4), 10) and 11), we call C(x) the hat-Lagrangian (density). One can 
choose any value for the time xq = 7^ n , when interaction-picture fields and Heisenberg- 
picture ones coincide. As in, e.g., Refs. 2) and 3), we choose T in = — oo, which we call 
the initial time. It is well known that the theory is ultraviolet (UV) renormalizable. As 
a matter of fact, in general, introduction of renormalization counter terms in vacuum 
theory is sufficient for the out-of-equilibrium theory to be renormalized. Thus, in the 
following, we do not write down explicitly the UV-renormalization counter terms. The 
free hat-Lagrangian Co, being bilinear in the fields, is obtained from Q2.2|) and (]2T 



C = -ft{d 2 + m 2 )f 3 l (2.6) 
where % is the third Pauli matrix and 
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Then, we are naturally lead to a two-component theory: The propagator, the vertex 
and the self-energy part take 2x2 matrix form. We use B to denote the 2x2 matrix 
whose component is Bij. From ( [2.5| ), we see that the vertex matrix V has simple 
structure, V = diag(t> , —v) with v the vertex factor in vacuum theory. The propagator 
fl2.3|) is translated into the propagator matrix: 



A(x,y) = -i{T c (j>(x)ft(y))), (2-7) 

where Tq acts as rearranging the fields in time-path order by recalling the above def- 
inition of the type-1 and type-2 fields. In this and the next sections, taking Co (Eq. 
(|2.6|) ) as the free hat-Lagrangian, we construct a perturbative scheme, which we call the 
bare-iV(umber) scheme. 

Fourier transforming the propagator A(x,y) [the self-energy part E(x,y)] on x — y, 
we have A(X; P) [t(X; P)} (cf. Eq. ([LI])). From the setup in §1, A(X; P) and S(X; P) 
vary slowly in X. Then, we employ the derivative expansion: 



G(X;P) 



1 + (X - Yfd Y » + ^{X - YY(X - Y) v d Y .d Yv + ... 



G(Y;P), (2.1 
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where G stands for A or E. Throughout this paper, unless otherwise stated, we keep 
only first two terms in (^.8|) (gradient approximation): 



G(X; P) ~ G(Y- P) + (X-Y)- d Y G(Y; P), 
where and in the following, '~' is used to denote the gradient approximation. 

2.2 Quasiparticle representation for the propagator 

In this subsection, pigeonholing the standard approach, [TJ, 01 H, ^ [TO] we introduce a 
quasiparticle representation for the propagator. 

2.2.1 Preliminary 



From (|2.6| ) and fl2.7|) , we have 

(d 2 x + m 2 )A(x, y) = (8 2 + m 2 )A(x, y) = -f 3 5\x - y) 

and from ( |2.7| ), we can readily obtainf] 

A n + A 22 = A 12 + A 21 , (2.9) 

A R = A n - A 12 = -i6{x - y ) [<j>(x), 0%)] , (2-10) 

A A = A n - A 21 = -i6(y - x ) [<j>(x), <P ] (y)] , (2.11) 

A c = A 12 + A 21 = -i (<f>{x)cf>\y) + ^(y)(f>{x)), (2.12) 

where Ar, Aa and A c are retarded, advanced and correlation functions, respectively. 
Equation Q2.9|) tells us that out of four elements Ay = 1,2), only three are inde- 
pendent. For independent quantities, we choose An, Aa and A c , which satisfy 

(d 2 x + m 2 )A R{A) (x,y) = (d 2 + m 2 )A R{A) (x, y) = ~5 4 (x - y), 

(d 2 x +m 2 )A c (x,y) = (d 2 + m 2 )A c (x,y) = 0. (2.13) 

It is to be noted that ( ^9|) is valid also for self-energy-part resummed propagators (see 
Appendix B). 

*Note that [4>{x), <ft{y)\ is a c-number function. 
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2.2.2 Quasiparticle representation 

Now we introduce a short-hand notation F ■ G, which is the function whose u (x,y) 
component" is 

[F-G](x,y) = [ d 4 zF(x,z)G(z,y). 



Motivated by equilibrium thermal field theory (ETFT),[Q, ^, [LI], [12] we introduce new 
functions [TJ, § [TO| 

/b ( x > y) = s(x - y ) g B fa y; x o) 

(with 'S' for short of Bare- A" scheme) and write A c (x,y) (= A c + Ac ) in the form, 

AW = A? • (1 + 2/W) - (1 + • A?, (2.14) 

where A( + )'s [A^'s] stand for the positive [negative] frequency part of A's (cf. ( |1 . 1| ) ) 
and '1' is the function whose u (x,y) component" is 5 4 (x — y). (The representation of 
A c in terms of a vacuum-theory kit is given in Appendix C.) As in ETFT, in the case 
where p is diagonal in momentum space, A(x, y) = A(x — y) and 

d 3 xe~^g^fa = rn T (p) - 9(-r) (r = ±) 



with h + the number density of the quasiparticle and n_ the anti-quasiparticle (cf. 
Appendix C). In the present out-of-equilibrium case, g^(x, y; x ) depends on X (= 
(x + y)/2) and on x only weakly. Then, in what follows, for g^(x, y; Xo), we use the 
gradient approximation not only for X (= (x + y)/2) but also for Xq. 

It is worth noting that A may be written as (cf. fl2.9|) - Q2.12| ) and (|2.14|) ) 



T = ± 

A d ia 9 = ^A^ = diag(A fi , -A A ), 



=± 

/ 1 1 A / f W\ _ ( 1 + f i B 



BL(fS') = , , B R (A r ')= . . • (2.16) 



This representation may be interpreted in terms of "retarded/advanced-quasiparticle 
picture." We introduce "retarded/advanced-quasiparticle fields," 4>ra and 4>rai through 
transformations (with obvious notation), 

* = E ^ (T) = E bmP) ■ & f = E ^ (T)t = E • ^(/b } )- (2.i7) 
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From (|2.16|) , we see that 



B L (f^).lT 3 . B R U %>) = lf 3 , 
which guarantees the canonical commutation relation to be preserved: 



0(x), (y) 



S(x - y ) 



■J 

•ra(x), <p RA (y) 
if 3 5 4 (x-y). 



5(xo - Vo) 



Thus, the ETFT counterpart of ( [2.17| ) is called the (thermal) Bogoliubov transformation. 
il| The statistical average of T c \'Pra4 > ra) ( c ^- (P-7P) assumes the form 



(T C (<t>RA(x) 4>R A (y))) = iA daig(x ~ y) . 



(2.18) 



Thus the fields 4>r A and (j) RA describe well-defined propagating modes in the system. It 
is to be noted that 4> RA is not the hermitian-conjugate of 4> RA , which is a characteristic 
feature of the theory of this type||, 0. Substituting ( gT7|) and ( g3§ into Q, QZg) - 
( |2.12j ) with (|2.14|) are reproduced. 

Incidentally, nonequilibrium thermo field dynamics [f|, [l(| is formulated by taking 
( p. 17] ) as one of the starting hypothesis. 

Now let us go to the P-space (cf. (|OJ)). A R and A A go to (cf. ( gHj ) and ( gTTp ) 

1 



A 



p2 _ m 2 _)_ ^ e (p Q )0+ 

It is straightforward to show that (|2.14| ) goes to 

A C (X;P) ~ -2me(p )[l + 2g ( i {po) \X;p)]5(P 2 -m 2 
+2 i d x »g« P0)) (X;p)d P ^ 

where 



(2.19) 



m 



2 • 



(2.20) 



9 B ±) (X;p) = j rf(x-y)e-P( x -^ ±) (x,y;X 



with X = (x + y)/2. One can easily show that (|2.13|) with Q2.14Q yields 

d x J±\X; p) = p ■ V x g£\X; p) = 0. 
Then, the last term in (|2.20|) vanishes:| 

A C (X; P) ~ -27TZ6(p )[l + 2^ (P0)) (X; p)] 5(P 2 - m 2 ). 



(2.21) 



(2.22) 



^As has been mentioned in §1, ( |2.22| ) is meaningful in the region \P»\> If Li 1 , where (|IJ) makes 



2.2.3 Particle-number density 

In order to find out the physical meaning^] of , let us compute the "free" or nonper- 
turbative part of the current density: 



fix) 



<P\x) & 1 4>{x)-4>{x) <9 M <p\x) 



where d= d— d- Taking the statistical average of these quantities, we obtain 



(f(x)) = -(d^-d x JA c (x,y) 



(2.23) 



Straightforward manipulation using Q2.22 ) yields 

d 3 p 



(f(x))-(0\f(x)\0) 
(j(*)> - (Ojj(a;) 0) 



(2vr) 3 
d 3 p 



g£\x;p) + {l + g ( B >(x; p)} 



(-)/ 



P 



{l + ^fop)} , (2.24) 



where E„ 



(2tt) 3 E p 

\/p 2 + m 2 . Equation Q2.24j ) is to be compared with 

d 3 p 



(f(x))~(0\f(x)\0) 



(27T) 



3 V B 



N^\x;p)-N ( ->(x;p)} 



(-)/ 



(2.25) 



where v B 



[l,p/E p ) is the four- velocity and N B (x;p) 



[N B \x;p)] is the number 



density (function) of the quasiparticle [anti-quasiparticle] with momentum p. Comparing 
(pi) and (ggp , we have 



N, 



(+), 



B 



x;p) 



9 B \ x ;p) 



N ( ~\x 



sense. (In the region |P**| < l/L^, A ( c ±] in ( |2.14|) "mixes" with A 



p). (2.26) 
For dealing with the region 



|_P M | < 1/£ M , one should return back to the fundamental form (2.12) or rather to (C.2) in Appendix 



C. Nevertheless, in most practical cases, one can use the form ( 2.22 ) for whole range of P^. This 
can be seen as follows. Let T be a typical scale(s) of the system under consideration. In the case 
of thermal-equilibrium system, T is the temperature of the system. Due to interactions, an effective 
mass M inc [(X) is induced. In the case of m >> \f\T, Mi n d(X) is not much different from m and, for 
to < VxT, a tadpole diagram induces M lnd {X) of O(VAT). [Since V\T ( or even AT) is the scale that 
characterizes microscopic reactions, the setup in §1 shows that l/L^ « \f\T (or AT).] Since most 
amplitudes, when computed in perturbation theory, are insensitive to the region \P^\ << Mi n d(X) of 
A C (X;P), the precise form of A C (X;P) in this region is irrelevant. In the following, we keep this in 
mind. Incidentally in the case of equilibrium thermal QED or QCD (to = 0), there are some quantities 
that diverge at infrared limits to leading order of hard-thermal- loop resummation scheme |l3|, [l2] |. For 
analyzing such quantities, one should use the fundamental formula (2.12) in the region |P^| < l/L^. 
^See also Appendix C. 
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It is to be noted that the argument x here is a macroscopic spacetime coordinates, which 
we write X throughout in the following. 



2.2.4 \po\ prescription 



Let us write fl4|| N B (X, ±p) = N B ^'(X; E p ,±p), where p = p/|p|. In the case of 



r(±) 



ETFT, N B ±] = N ( b>(E p ). From (^22|) , wc sec that gf\X, p) appears in combination 



r(±) 



(±)/ 



with 5(P 2 



5(Po — Ep). Then, at first sight, it seems that no difference arises 



between N^(X; E p , ±p) and N B ± \X; ±p , ±p). It is well known in ETFT that this 
is not the case, since, in general, A C (P) is to be multiplied by the functions that are 
singular at |po| = E p . The correct choice has been known [|15f| for some time now to be 



N B ± \±p ) (= N B ^'(\p \) — the "|po| prescription." On the basis of this observation, we 
assume^ that the |p | prescription should be adopted: 

g ( B +) (X;p) = N B + \X;E p ,p)^N B + \X;p ,p) = f B + \X;P), 



r(±) 



9 ( b\X;p) 
Ib(X;P) 



1 - N B -\X; E p , -p) - -1 - N [ B '(X; -p , -p) = fe>{X; P), 



<■(-), 



9(p ) f B + \X; P) + 9(-p ) f B ->(X; P). 



(-)/ 



(2.27) 



It is to be noted that the translation into the |po| prescription is formally achieved 
by (cf. - fl2Tp ), 



A = 

where fs(x, ?/) is defined by 



A h -(1 + 2/b)-(1 + 2/b).Aa, 

BlUb) ■ ^diag ■ B R (f B ), 



(2.28) 



fii{x,y) 



d 4 P 

WT 4 



-iP-(x-y) 



f B (X;P) {X = {x + y)/2) 



with f B (X; P) as in ( gjg) . Now (jjTgjJ) turns out to 

P-d*/ B (X;P) = 0. 



(2.29) 



This is a continuity equation for along the "flow line" in a four- dimensional X M -space. 
Equation (2.29) may be solved as 



f B {X- P) = f B (X = T m , X - (X - T in )p/ Po ; P), 



(2.30) 



"it can be shown that, as in ETFT, ( |2.27| ) is consistent with mass- derivative formulajl^, [l5| at least 
up to the terms with second-order X^-derivative. 
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where T in (= — oo) is the initial time (cf. §2.1). The propagator A(X;P) with this 
f B (X; P) takes the form (cf. (ggg)), 

A C (X; P) = -27TZ e(po)[l + 2/ fl (X; P)] 5(P 2 - m 2 ). (2.31) 

It is to be noted that ( |2.29| ) guarantees that all higher-order terms in the derivative 
expansion (cf. ( |2.20|) ) vanish and, in this sense, (|2.31|) is exact. 



3 Bare- TV scheme 

Interactions among the fields give rise to reactions taking place in a system, which, in 
turn, causes a nontrivial change in the number density of quasiparticles. In this section, 
within the bare-iV scheme, we analyze the self-energy-part resummed propagator and 
derive a generalized Boltzmann equation. 

3.1 Self-energy part 



We start with noticing that Y^i - =1 Ey (x, y) — 0, a proof of which is given in Appendix 
B. Using this property and (|2.16|) , we obtain 



B R {f B ) ■ E • B L (f B ) = ( S * ^ ) , (3.32) 

E_r = En + E12, Sa = — (S 2 2 + Ei 2 ), (3.33) 
^off — E12 + E12 • f B — E 2 i • f B + E^ • f B — f B ■ Ea, (3.34) 

where E# and E^ are called the retarded and advanced self-energy parts, respectively. 
Going to momentum space, we obtain 

E o// (X; P) ~ i {f B (X; P), PeE fi (X; P)} + if P), (3.35) 



where 



(A( Y .P^ R(Y D „- dA(X;P)dB{X;P) dA(X; P) dB{X; P) 
{A(X,P),S(X,P)} = — — — — , (3.36) 



1 



f P) = (1 + / B (X; P))E 12 (X; P) - / B (X; P)E 21 (X; P). (3.37) 



In deriving ( |3.35| ), use has been made of the relation E j4 (X; P) = [E B (Jf; P)]*, a proof 



of which is given in Appendix B. 
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3.2 Self-energy-part resummed propagator 

Let us compute the self-energy-part resummed propagator, with the aid of the Schwin- 
ger-Dyson equation, 

G = A + A-£-G = G-£-A + A. (3.38) 
Using flggj ) and (|3T32D , we obtain 

G = B L (h) ■ ( G * ~° R ^ ' ° A \ ■ B R (M, (3.39) 

Gr(A) = ^R(A) ■ [l + ^R(A) ■ Gr{A)] ■ (3.40) 

As a matter of course, G's obtained from ( 3.39 ) are written in the form (|2.9|) - (|2.12 ) 
with G R (G a ) for A R (A a ) and, for A c , 

G c = G R ■ (1 + 2/ B ) - (1 + 2/ B ) ■ G A - 2G R ■ E o// • G A . (3.41) 

Let us compute G(X;P)'s (X = (x + y)/2), the Fourier transform of G(x, ?/)'s oni-t/, 
to the gradient approximation. In ETFT, ^ ff = 0. Furthermore, as will be shown 
below, ^ ff is proportional to dN^ /dX^, where are (the contribution to) the 

physical number densities. Then, in the derivative expansion for T, ff(x,y) (cf. (|2.8|) ), 
we keep only zeroth order (no-derivative) term. Thus, we obtain 

1 

P 2 _ m 2 _ E R{A JX; P) ± ie(p o )0' 
G C (X;P) ~ G?\X;P) + G^(X;P) + G%(X;P) : (3.43) 

where 

G^(X;P) = [l + 2f B (X;P)][G R (X;P)-G A (X;P)], (3.44) 
G$(X;P) = -{f B (X;P),ImX R (X;P)}(G 2 R (X;P)-G 2 A (X;P)) 

-if P) P) + G 2 A (X; P)) , (3.45) 

GS(X;P) = E off (X;P)(G R (X;P)-G A (X;P)f. (3.46) 

It is to be noted that ( |3.42| ) is exact to the gradient approximation, i.e., no term including 
first derivative (with respect to X^) arises. 

In narrow-width approximation, ImY>R = -/mS^ — > —e(p )0 + , GrGa involved in 
G^ develops pinch singularity in a complex p -plane, while G^ and G^ turn out to 
the well-defined distributions. 



G R{A) (X;P) ~ - — 7 -— r — 7 - T - z , (3.42) 
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For clarifying the physical meaning of G^ 2 > we compute the contribution to the 
physical-number density through analyzing the contribution to the statistical average 
of the current density (f), (|2T23| ) with G c for A c . G^ as well as G ( c ] - A c lead to 
perturbative corrections to in (|2.24j ) due to quantum and medium effects, while 
G^*2 yields the large contribution, 

d 4 P 



(5f(X)) = i 



(27T)' 



P» E o// (X; P) (G R (X; P) - G A (X; P))' 



(3.47) 



which diverges in the narrow-width approximation. In the next subsection, we inspect 
this large contribution more closely. 



3.3 Generalized Boltzmann equation 

The contribution of G^ to the physical number density, SN^\ is obtained from (|3.47|) 



(cf. Q2.25 )), which we carry out in Appendix Dl: 



57V (±) (X;±cj ± ,±p) 



2 llmE^I 



PO 



± 



2u± 



[1 + N^)Z 12{21) 



1 



2co 



{f B , ReZn} 



± 



'£21(12) 



J p =±U)± 



(3.48) 



(3.49) 



Here u± (= lu±(X; ±p)), Eq. (|D.1| ) in Appendix D, is the energy of the ± (quasiparti- 
cle/anti-quasiparticle) mode with momentum ±p, Z± (Eq. ( p.3| )) is the wave-function 
renormalization factor and t± = 2o;±/|/mEfl(X; ±o>-t-,p)| is the time during which the 
± mode damps. T±\ Eq. ( |3.49| ), that has come from G — «£ // in ( |3.37| ), is the 
net production rate of the quasiparticle (anti-quasiparticle) of momentum p (— p). In 



fact, r 



(p) 



(p) 



is the difference between the production rate and the decay rate, so that r ± 



is the net production rate. As mentioned above, in the case of equilibrium system, 
T±* = (detailed balance formula). Note that ( |3.48j ) diverges in the narrow- width 
limit, t± — > oo. 

Recalling the notion of the spacetime cells (cf. §1), in Appendix D2, we derive from 
( [3. 48] ) a generalized Boltzmann equation for = N^\X; ±uj±, ±p) + 
5N^(X; ±o;±, ±p) with Ng the bare number density in §2: 

r(±), 



dX n 



+ v± • V x iV (±) 



Z ± Y [ l ] 



du ± dN^(X;±p ,±p) 



dX tl 



dP» 



(3.50) 



p =±UJ± 
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The first term on the right-hand side (RHS) is the collision term, while the last term 
represents the effect due to a change of "mass." Discussion of ( |3.50| ) will be postponed 
until §4. 

The physical-number-density J\f^ may be written as 



Here, as mentioned above, AN^\ being the contribution from G c +G : 
is a "regular correction." 



r cl 



(3.51) 
A, (e GX 



3.4 Absence of large contributions 

In perturbatively computing some quantity, G C (X; P) appears in the form, 

d A P 



(27T)' 



G C (X;P)F(P), 



(3.52) 



where F{P) is the function that is reg ular at p = ±u ± . G { c ] + G^ G G c , Eqs. (|3~43p - 
( p.45| ), leads to a "regular contribution" to fl3.52|) . Equation ( |3.52|) with G^ + for 

G c reads 

d 4 P 



G i ?\X-P) + G^(X-P) HP) 



which diverges in the narrow-width limit. Using ( |D.4| ) (in Appendix D) in G 
(G R - Ga) 2 S o// ), we obtain 



(3.53) 

(i) / 

c2 l — 



Eq. (3.53) 



E 

r=± 
+ .. 

E 

r=± 
+ .. 



d 3 p Z± 
(2tt) 3 2uj± 

d 3 p Z± 
(2tt) 3 2uj± 



±[l + 2f B ± \X;±u ± ,p)} 



zE 



(±) 

off 



1/mE 



(±)i 

R I 



1 + 2 \zi l N {±) + (1 - Zj^iV 



(±) 



(3.54) 
.F(±w±,p) 

(3.55) 



where stands for "regular contributions." In obtaining (|3.54|) , we have assumed that 
the quasiparticles here are well defined, i.e., JmS^ 



(±) 



7mE^|. Thus, in terms of 



there does not appear the contribution that diverges in the narrow- width limit. 
The "regular contributions" (|3.55|) are a functional of Jb, which, if one wants, can be 
rewritten in terms of or A/"^ by using the solution to (|3.50|) and (|3.51|) . Compu- 
tational procedure, within the bare-A" scheme, of a generic amplitude is similar to the 
one within the physical- A" scheme to be explained in §4. 
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4 Physical- N scheme 



4.1 Preliminary 



Introducing new functions f^(x, y) = 5(x$ — yo) (/^(x, y; xq), we redefine the fields, (/> 
and 0t, by ( gTffl with B L (r)U {t) ) B L{R) {f^ ] ) (cf. (gl6|)). Then, the fields in this 
scheme are different from the fields in the previous sections. In fact, that / does not 
satisfies ( [2.21 ) or ( |2.29| ) means that ( |2.13[ ) does not hold. Then, the interaction-picture 
field <p and then also $ do not obey the Klein-Gordon equation, (d 2 + m 2 )<j) ^ 0, which 
means that Cq in (|2.6|) is not the free hat-Lagrangian. For the purpose of finding the 
correct free hat-Lagrangian C' , we proceed as follows: Compute (d 2 + m 2 )A(x,y) with 
A(x,y) as in ( |2.15| ) with Q2.16Q (with —* /^)- The result may be written in the 
form, 



(r) 



T=± 



D T {x,z) = f 3 {d 2 x + m 2 )5\x-z) + A[g< T \x,z; 



x 



d 



+2<jM(x, z; x )— - (V 2 - V 2 z )g {T \^ z; x ) 

OXq 



6(x - Zq), 



where g(x, z;x ) = <9g(x, z; x )/dx , etc., and A is a matrix with A^ = (— 
1,2). This means that the free hat-Lagrangian jC' q is 



4 = - / rf 4 ^0 (r)t (^)^(^y)0 (r) (y)- 



(4.56) 



Starting with (<9 2 + m 2 )A(x, y), we are led to the same C' as ( |4.56| ) above, as it should 
be. Since C ^ Co, Eqs. ( |4.56|) and ( |2.6| ), the counter term C c = C — C f appears in the 
hat-Lagrangian, which yields the two-point vertex function, 

d 



iV(x,y) 



T = ± 



^(x,y;xo) + 2^(x,y;x 



dxo 



-(V 2 x -V 2 y )g( T \x,y;x )] 5(x -y ) 



d 4 P 



-iP-(x-y) 



^(rpo) [P-d x g^\X;p)} A, 



T=± 



where X = (x + y)/2. Going to the |po| prescription, we have (cf. ( |2.27| )' 



iV(x,y) (= iV{x,y)A) 



d A P 



AP-(x-y) 



[P-d x f{X-P)]A. 



(4.57) 
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It should be emphasized that this result is exact, not relying on the gradient approxi- 
mation. 

The propagator A(X; P) takes the form as in §2, provided that A C (X; P) is given in 
( p.20| ) with — > Going to the |po| prescription, it reads 

A C (X;P) ~ -2me(p )[l + 2f(X;P)}5(P 2 - m 2 ) 
-2iP-d x f(X;P) [A* (P) + Aj(P)] 

( = Af+AP). (4.58) 

Now, E = B R ■ E ■ B L takes the form fiT5% ) with 

E o// = E 12 + E 12 - / - E 21 • / + X A ■ f - f - E A + V. (4.59) 

Fourier transformation leads to 

iKff(X; P) ~ {/(X; P), P 2 - m 2 - PeE^X; P)} - f»(X; P). (4.60) 



The E-resummed propagator is written in the form (|3.42 ) - ( p. 46] ) with /b — ► /. 
So far / has not been specified. 

4.2 "Renormalization condition" and the generalized Boltz- 
mann equation 

We now define /, such that the number densities 

NM(X;u + ,P) = f(X; Po = u; + ,p), 
N<>-\X',u-,--p) = -l-/(Jf;po = -w+,p), (4.61) 

are as close as the physical- number densities Af^. As seen in §3, i n ( |3.46|) , which is 
proportional to E jj, yields a large contribution to the physical- number density. Then, 
as the determining equation for /, we adopt E yj(x,y) = or 

E o// (X;P)=0, (4.62) 

which we refer to as the "renormalization condition" for the number density. It should 
be emphasized here that this condition is by no means unique. E G // is a "functional" 
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B R (f)-t-B L (f) 



of /, £<,//[/]• We can adopt £ //[/] = (/ ^ /), provided that f(X;p = ±u±,p) = 
f(X;p = ±w±,p). 

In the present scheme with ( |4.62| ), transformed self-energy part £ (cf. Eq. ( p.32|) ) , is 
diagonal: 

' Z R 

o -s A 

which is in conformity with the quasiparticle picture. Thus, the present scheme is similar 
in structure to ETFT: Fourier transform of A • £ • A, etc., are free from pinch singularity 
and then no pinch singularity appears in perturbative calculations. Similarly, G C (X; P) 
and then also G(X; P) are free from pinch singularity in the narrow-width limit (cf. 
above after Q3.46Q ). 

Let us inspect the physical implication of the "renormalization condition" ( 4.62| ). We 
write ReE R (X;P) = M 2 (X) + ReE' R (X; P), where M 2 (X) is the contribution from a 
tadpole diagram. Note that the tadpole diagram does not contribute to T^ p \ Equation 
( OT) with becomes 



v ' ; dX» dP^ 
df(X; P) d[M 2 (X) + ReE' R (X; P)} 



dPf, dX» 



(4.63) 



Recalling that JV<+> = JV«(X;po = u;+,p) [N^ = N^(X; -p = cj_, -p)] in (TO) 
is the number density of the quasiparticle [anti-quasiparticle] with momentum p [— p], 
it is straightforward to show that ( |4.63| ) becomes, on the mass-shell, po = ±u;±, 

1 v ± • V X A^ (±) ± 



I po = ±CJ± 

dN^(X;uj±(X;±p),±p) du± dN& du± dN& 



dXo dp dX. <9X dp 

(4.64) 

po = ±w± 



z+r (p) 



which is (almost) identical to ( |3.50| ) in the bare-iV scheme. iV^) = N^(X; u±, ±p) 
here is essentially (the main portion of) the relativistic Wigner function and ( |4.64j ) is 
the generalized kinetic or Boltzmann equation (cf. Ref. 17)). 
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4.3 Procedure of solving ( |4.t)3| ) 



Equation (|4.63|) determines f(X; P) self-consistently. One can also approximately solve 
( 4.63Q , order by order, through the following iterative procedure. At the RHS of ( [4.63 ), 
three small quantities^ are involved; the coupling constant A, r( p ) and df /dX^. T^ pS) 
measures how the system is far from equilibrium. P^dx^f has been determined by 
( |4.63| ) at lower orders (than the order under consideration) but the other components of 



dxt^f reflect the form of initial /| 



B 



For simplifying the following 
presentation, we assumep] that the interplay of all these small quantities yields three 
terms, on the RHS of ( f4.63| ), which are of the same order of magnitude. [For other cases, 
(trivial) modification of the following procedure is necessary] Let us write the RHS of 
( |4.63| ) TZ[A c (f)}, which is a functional of /, and expand /, A c and 1Z (with respect to 
the small quantity): / = £\ =0 /«. A c = £\ =0 A®, TZ = £\ =1 K®. Here Af are as in 
(g3|) and /(°) is nothing but f B in §2 - §3 and satisfies flggg ), P ■ d x f {0) = 0, whose 
solution is given in fl2.30p . 

Zeroth order: Equation (|463|) reduces to (|2~29|) : P ■ d x f {0) = 0. 
First order: Equation ( |4.63| ) becomes 



IP ■ d x f {1) (X; P) = ^ (1) [A( 0) (/ (0) )]. 



(4.65) 



Solve ( |4.65| ) under the condition | = ^ =0. 

Second order: Compute T, ff(X\P) up to second X M -derivative terms in the derivative 



expansion, which yields the additional term on the RHS of ( |4.63j ). 



Then, we have 



«9 2 (£ 



12 



J 21, 



d 2 f <9 2 (£ 



dX v 

<9 2 (£l2 — S 2 l / 



dP^ dP v 
d 2 f 



+ 



12 



J 21, 



d 2 f 



OP^ dP u 



dX^ dX" 



dX^ dP v 



dP^ dX» 



P-dxf 



(2) 



7e(3)[A(o)(/W)] + {ftW[A(°)(/(°) + /(D) + A«(/(°))] 
-KW[A®(f<®)]} +K'[AW(fW)\, (4.66) 



"Note that dReY, R /dX is a functional of df/dX. 

"For the case of m 2 < AT 2 (with T the typical scale(s) of the system) the tadpole self-energy-part 
resummed propagator^, ^] should be substituted for the bare propagator A. 
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Since P ■ d x f {0) = 0, A ( c ] [f {0) ] = (cf. flOSD ). Solve flE5E|) under the condition 

Higher orders: Proceed similarly as above. Further "improvements" of ( |4.58| ) and (|4.63|) 
are necessary. 

More efficient way is to construct S-resummed propagator, Eqs. ( |3.42|) - (|3.45 ) 
(So// = 0), at each order and use it for yet higher-order calculations (cf. below). 

4.4 Procedure of perturbative calculation 

Perturbative computation, to iVth order, of an n-point amplitude A{P\, ...,P n ), 

/n 
[Jrf 4 a ; i e- iE ' P, ' I %,..,^) ) (4.67) 
i=i 

goes as follows. 

1) Draw all relevant skeleton diagrams. Contribution to A of each skeleton diagram 
may be written in the form, 



/ Y[d A ZiS(x l ,...,x n ;z l ,...,z j ). 



(4.68) 



Here the external-vertex (spacetime) points and the internal- 

vertex (spacetime) points, which we collectively write £i, ...,£ n+ j. 

2) For the propagators involved in each skeleton diagram, use the self-energy-part 
(S) resummed propagator G, with S and / computed to appropriate orders. [For the 
skeleton diagram that is already of iVth order, it is sufficient to use the bare propagators 
and the "bare" /'s.] 

3) To analyze £(£1, £n+i)> as (not necessarily independent) variables, employ the 
center-of-mass coordinates of A, X (= X^=i x */ n )> an< ^ relative coordinates £ ifc — £j k 
(k = 1, ...,m with m a number of propagators). Here £j fc and £j k in ^ k — ^ are the 
points that are connected by the propagator (j(£i fc , £j fc ) (= Gt) , 

G k = G (^+^5 - = J ^ e-^<^G (£*±&; Q k "j , 

whose form is known above. [An external as well as an internal point appears, in 
general, in several relative coordinates.] It can easily be shown that (& fc + £j fc )/2 may 
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be "decomposed" as 



(4.69) 



where c's are numerical coefficients. The decomposition ( f4.69| ) is not unique, which is 
not the matter for our purpose. Carry out the derivative expansion: 



d ® k e -iQv{Zi h -ti k ) 
(2vr) 4 



1+ J2 c ^-^- dx+ - 



G{X;Q k )). (4.70) 



Pick up the term with — £ J; ). Multiplication of Gi (e S) yields 



c k i^ l -^ jl )-d x G(X;Q k )G l 



c H d x G(X;Q h 



d A Q 



1 d 



G 



ji . 



Qi 



(27f) 4 •■ / ! ••• _ 

The higher X M -derivative terms in ( 4.70|) may be dealt with similarly. Deal with 
all the propagators in the above manner. For /(X;Qz)'s, substitute the solution (of 
appropriate order) to (|4.62|) . 

4) Carry out the integrations over #'s in (|4.68 ), which yields j momentum-conserva- 
tion ^-functions at the internal points, z's. 

5) Carry out the integrations over j Q's. Add the contributions of all the diagrams 
to obtain A(x\, ...,x n ), which, as is obvious from the above procedure, is of the form 
A(X; xi - x 2 , x n _i - x n ). 

6) Substitute A thus obtained into (|4.67|) to obtain 



A(P x ,...,P n ) = j d 4 XA(X;P 1 ,...,P 1 

A(X;P u ...,P n ) = ! 

J i=i 

x A(X] xi - x 2 , x n _i - x n ) 



(4.71) 



H d\x l+l - x^ e-^i PH** + i-*i) 



(4.72) 



where P[ is the linear combination of Pi, P n . Carry out the integration over — x^ 
(i = l,...,n — 1) to obtain A(X; P\, P n ). Note that the momentum conservation 
Yli=i Pi = holds (cf., however, the footnote in §2.2.2). Note also that A depends 
weakly on X through /(X; Q/J's. Were it not for this X-dependence, integration over 
X in ( pTTTD would yields (2tt) 4 5 4 (^ PA. 
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5 Comparison with related work 



As has been mentioned in §1, the physical- N scheme is employed in the literature. 
5.1 CTP formalism 

The derivation of the generalized Boltzmann equation (GBE) in Refs. 2), 7), 18) and 
19)0 starts with the Schwinger- Dyson equation, Eq. ( |3.38D , 

(A" 1 - e) G = G (A" 1 - e) = 1. (5.73) 

Here we stress that this equation is nothing more than the equation that serves as 
resumming the self-energy part to makeup the resummed propagator, (|3.40|) and ( |3.41|) . 



Then in order to derive the GBE, an additional input or condition is necessary. Our 
condition is the "renormalization condition" £ // in ( |4.62j ). 

We first see the derivation in Refs. 7), 18) and 19), in which the lowest nontrivial 
order is dealt with. The condition adopted there is essentially 

G C (X; P) = [l + 2f(X; P)} [G R (X; P) - G A (X; P)}. (5.74) 

[To distinguish from ours, we write /' for / here.] No counter Lagrangian is introduced, 
so that the solution for G c to Q5.73j ) is given by ( |3.41| ) with f B — > f (cf. the argument 
in §4): 

G c = G R ■ (1 + 2f) - (1 + 2/') ■ G A - 2G R ■ E o// ■ G A , (5.75) 

where E D // is as in ( |3.34|) with f R — > f ■ Note that the condition Q5.74! ) is the leading 
term of the (Fourier transform of the) exact form ( |5.75 ) with E j/ = o. 



Fourier transforming ( |5 . 75[ ) and equating with ( |5 . 74| ) , one obtains 
{f'(X; P), G R (X; P) + G A (X; P)} ~ -2iG R (X; P)E off (X; P)G A (X; P) . (5.76) 



To the lowest nontrivial order under consideration, the left-hand side (LHS) of ( 5.76Q is 
approximated as 

LHS of (5.76) ~ -2iP-d x f'(X;P){G R (X;P) + G 2 A (X;P)} 

= AiP ■ d x f(X;P)G R (X;P)G A (X;P) ( Po = ±cu ± ) . (5.77) 



ttln Ref. 19), nonrelativistic many-body theory is dealt with. 
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Using ( |5.77| ) in ( |5.76| ) on the mass-shell, po = ±u±, one gets the Boltzmann equation. 



In Ref. 2) is given an all-order derivation of the GBE without introducing the counter 
Lagrangian. Although the correct GBE comes out, the derivation is not self consistent. 
In order to see this, we first note that, upon using (|3.33| ), (|3.34j) and (|3.40|) , the solution 
( |5.75| ) may be rewritten as 

G rG c G = G\ + (j2, (5.78) 
G x = Eu + Eaa (5.79) 
G 2 = 2(f-A A 1 -A R 1 -f). (5.80) 

At an intermediate step of solving Schwinger-Dyson equation, there appears the term 
proportional to A^ 1 • A c • A^ 1 , which turns out to G 2 - In Ref. 2), G 2 is missing, and 
thus, when interaction is switched off, G c does not reduce to the free form ( p.28| ). [It 
is to be noted that, in ETFT, Fourier transform of A^ 1 ■ A c • A^ 1 is proportional to 
(P 2 — m 2 ) 2 5(p 2 — m 2 ) and then vanishes.] As the additional condition, in place of (|5.74|) , 
G c = Gr-(1 + f) — (l+f')-Gj[ is adopted. This again contradicts with the exact solution 
fl5.75j ) unless E // = 0. By equating two relations above, the correct GBE comes about. 



As has been shown in §4, for constructing a consistent physical- N scheme, the counter 
Lagrangian C c should be introduced. As mentioned above, however, C c is not introduced 
in Refs. 2), 7), 18) and 19), so that no self-consistent perturbative calculational procedure 
is presented. 

5.2 Thermo field dynamics (TFD) 

In contrast to the CTP formalism, the TFD counterparts of <pi and <p2 are independent 
fields and mutually commutable. 

We first consider the equilibrium case. In spite of the fact mentioned above, as far as 
the perturbation theory is concerned, TFD is equivalent to CTP formalism. However, 
we should mention one point — the |po| prescription (§2.2.4). In the CTF formalism, the 
|po| prescription is deduced [13] from first principles. The very existence of the vertical 
segments^], 12, 15 1 of the time-path in a complex time plane plays an important role. In 
TFD, however, there is no counterpart of the vertical segments, and then, on the basis 
of the consistency argument (with the spectral representation and KMS condition), the 
\p \ prescription is assumed. [O, [L6| 
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The nonequilibrium TFDQ, (NETFD) is formulated on the basis of (nonequilib- 
rium generalization of) the so-called thermal-state condition. As a consequence of this, 
initial correlations (cf. ( |2.4|) and Appendix A) are absent. Thus, when compared to the 
(CTP) formalism developed in this paper, NETFD may be applied to more restricted 
class of nonequilibrium systems, i.e., the systems which are described by suitably de- 
fined quasiparticles, in terms of which the generalized thermal-state condition holds. On 
the other hand, NETFD is a "single- (space)time" canonical formalism, without distinc- 
tion between the microscopic- and macroscopic-(space)times. In this sense, NETFD has 
much wider applicability than the CTP formalism presented in this paper, the latter of 
which applies only for "out-of-equilibrium systems." Incidentally, in NETFD, there is 
one parameter (s-parameter), which has no counter part in the CTP formalism. 

Our condition for determining the "physical-number density," N^, is that, with 
pinch singularities (in the narrow- width limit) disappear. On the other hand, 
NETFD, in which the (space)time representation is employed, imposes an "on-shell 
renormalization condition," which results in the GBE. Since "the pinch singularity" is a 
notion in momentum space, it seems not to be immediately obvious how to translate this 
condition to the (space)time representation as adopted in NETFD. Nevertheless, closer 
inspection of the structure of both formalisms tells us that our condition is in accord 
with the "on-shell renormalization condition" in NETFD. Incidentally, how to reconcile 
the NETFD with the \po\ prescription, a notion in momentum space, remains to be an 
open question. 

6 Summary and discussion 

In this paper, we have dealt with perturbative framework of out-of-equilibrium rela- 
tivistic complex-scalar-field systems. We have assumed the existence of two different 
spacetime scales, the microscopic and macroscopic. The first small scale characterizes 
the microscopic correlations and the second large scale is inherent in the relaxation 
phenomena. 

We have proposed mutually equivalent two computational schemes, the bare-iV and 
physical-iV schemes. Both of them lead to the generalized relativistic kinetic or Boltz- 
mann equation. Then, we have presented the procedure of perturbatively computing a 
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generic amplitude. 

It is worth pointing out a similarity between two computational schemes presented 
in this paper and those in UV-renormalization in vacuum theory. Taking mass renor- 
malization as an example, let us see this. 

A "bare scheme" in vacuum theory takes Cq = —cf)\d 2 + m 2 B )<j), with tub the bare 
mass, as the free Lagrangian. Then, the propagator reads, A(P) = l/(-P 2 — m 2 B + i0 + ). 
Perturbative computation of the renormalized mass m yields m = Minis), which may 
not necessarily be the physical mass M.. Perturbative computation of some quantity, 
e.g., the physical mass A4, yields the result as a function of tub, which may be rewritten, 
using m = M(tob), in terms of m. The correspondence of the "bare scheme" here to the 
bare- AT scheme in §2 - §3 is as follows: C -> C in A(P) -> A(X; P) (Eqs. ( |2TTD|) - 
( CT ) with ( gig) and (gjjj])), m B -> JV^ , m -> iV^, m = M(m fl ) -> £7g. (gl8|) , 

A "physical scheme" in vacuum theory takes C = — 0^(<9 2 + m 2 )0, with m the renor- 
malized mass, as the free Lagrangian and then there emerges the counter-Lagrangian, 
£c = ^(m 2 — m 2 B )(j). The propagator is A(P) = 1/(P 2 — m 2 + i0 + ). m 2 — m 2 B is 
determined so that the perturbatively computed mass, M(m,m 2 — m 2 B ), is equal to m, 
m = M(m, m 2 — m 2 B ). The correspondence of the "physical scheme" here to the physical- 
N scheme in §4 is as follows: C -> C' in ( ^56|) , A(P) -> A(X; P) (Eqs. (P1D|) - (Pl2]) 
with (fTrj) and Q435|) ), £ c -> £ c = £ -A> m = M(m,m 2 -m|) -> E o// = in ( ggg ), 

To summarize, the bare-iV scheme is constructed in terms of the original bare-number 
density fs{X;P) [respect, the bare mass m B }. On the other hand, the physical-iV 
scheme is constructed in terms of the "renormalized" -number density f(X; P) [respect, 
the renormalized mass m]. Both perturbative schemes are equivalent. The first scheme 
starts with the "bare quantity" and the "renormalization" is done at the end, while in 
the second scheme, the "renormalization" is done at the beginning by introducing the 
counter hat-Lagrangian C c [respect. C c ]. It is worth recalling that the renormalized 
mass m is defined so as to absorb UV divergences. However, there is arbitrariness in the 
definition of the "finite part" of m, which is determined by imposing some condition. 
Under this condition, m is determined order by order in perturbation series. This is 
also the case in the present physical- N scheme. / (or N^) is defined so as to absorb 
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large contributions, which diverge in narrow-width limit. As pointed out above after 
( |4.62j ), there is arbitrariness in fixing the "finite part" of /. In §4, we have imposed the 
condition for determining /, under which f(X; P) turns out to be determined order by 
order in perturbation series. 



Acknowledgments 

This work was supported in part by the Grant-in- Aide for Scientific Research ((A)(1) 
(No. 08304024)) of the Ministry of Education, Science and Culture of Japan. 



Appendix A: Initial correlations 

Here we discuss to what extent, the initial correlations Cm-, Eq. (|2.4j) , may be ignored. 

We first note that, in the two-component representation (cf. §2.1), C n has 2 n compo- 
nents, all of which are identical. |2| Then, from Q2.10 ) and ( |2.11| ), we see that C2, which 



constitutes the medium part of the propagator A, does not appear in the retarded and 
advanced Green functions. This means that C2 does not appear in the two-point correla- 
tion in the linear response theory. With the aid of the standard formulae (cf. Sections 2 
and 5 of Ref. 2)), one can also show that C 2n does not appear in multi-point correlations 
in the nonlinear as well as the linear response theory. 
We decompose <fi(x) as in vacuum theory: 



/A j[ "** t " ,it * 1 ' (A - i) 



where qo = E q = \J q 2 + m 2 . a(q) [V(q)] in ( |A.1| ) is the annihilation [creation] operator 
of a particle [antiparticle] in vacuum theory: [a(q), a^(q')] = [6(q), 6^(q')] = 8(q — q'). 

We restrict our concern to C4. The argument below may be generalized to the case 
of general C2 n (n > 3). We are concerned only with the particle part of 0, the part with 
a(q) in ( |A.1| ), and of <p\ the hermitian conjugate of cf). Other parts may be dealt with 
similarly and the same conclusion will be obtained. 

Substituting ( |A.1|) and its hermitian conjugate to C 4 in ( p.4|) and carrying out the 
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Fourier transformation, we obtain, up to an numerical factor, 



2 



C 4 cx — 2 l — — - [*(p° - E Pj ) 5(q° - (q 1; q 2 ; Pl , p 2 ) c (A.2) 

llj=l V^Pj^Qj j=l 

(qi,q2;pi,P2) c = (qi,q2;pi,P2> - <qi;pi)(q2;P2) - (qi;P2>(q2;pi), 

where (qi, q 2 ; pi, p 2 ) = (a t (qi)a t (q 2 )a(pi)a(p 2 )) and (q; p) = (a f (q)a(p)), etc. The 
setup in §1, C 4 does not change appreciably in a single spacetime cell, leads to approxi- 
mate momentum-conservation: 

2 

c 4 ~ for - gy) > (a.3) 



Thus, we may write 



C, = A 4 ^(P, - g,) J] [5 + (Pf - m 2 ) 5 + (Q 2 - m 2 )] ,4, (A.4) 
\j=i J j=i 

where A(P M ) is the function whose width [height] is of 0(1/ L^ 1 ) [0(L M )] and satisfies 
/ A(P^)dPv = 1, and 5 + {P 2 - m 2 ) = 6(p )5(P 2 - m 2 ). 

A transition probability or a rate of a microscopic reaction is related to an on-shell 
amplitude. (In the case of equilibrium thermal field theories, this relation is settled in 
Ref. 20).) Let us now analyze the on-shell (P 2 = m 2 , etc.) amplitudes. As can be seen 
from ( |A.2| ), C 4 per se does not contribute to the on-shell amplitudes. Then, C 4 or C 4 
appears as a part(s) of an on-shell amplitude. Let us see the structure of a connected 
amplitude S, which includes a single C 4 . We may write, with obvious notation, 

S = [l[[d^d i Q,]C 4 T(Q 1 ,Q 2 ;P 1 ,P 2 ) 

d 3 pid 3 p 2 d 3 qi 

xAF(P u P 2 ; Qi, Q2 = Pi + P2 - Qi). (A.5) 



^ / n n n n A (£ pi + -^P2 ~ ~ ^|pi+P2-qi|) 



C4 is to be compared with the four-point Green functions G 4 's. For the sake of com- 
parison, we take the nth-order contribution, whose four legs are bare propagator 
with on-shell 5-function: 

oc A" +1 5 4 (j^{P 3 - QjU H [6{Ff - m 2 ) 5(Q 2 - m 2 )] B, 
\j=i J j=i 



26 



The counterpart of S, Eq. ( |A.5| ), is 



<j — A d(£/ pi + £, P2 - £ ?1 - ^| Pl+P2 _ qi |J 

J -C/pi rj P2 rj qi - C/ |pi+P2-Qi 

xBT(P 1 ,P 2 ;Q 1 ,Q 2 = Pi + P 2 -Qi). (A.6) 



From ( A. 5 ) and ( |A.6|) , we obtain, with obvious notation, 

5 (.4.F) 



(A.7) 



Noticing that ^4 and B are dimensionless, we assume that A and £> are of 0(1). Equation 
( |OD tells us that, unless {AT) / (BT) « 1/A n+1 , C 4 may not be ignored. 
Equation ( |C.2| ) in Appendix C, together with the setup in §1, tells us that 

(at(q)a(p)) ~ for |q-p|>l/|L|. (A.8) 

Motivated with this relation, we consider a case, where (qi, q 2 ; Pi, p 2 ) c has the property 
of factorizability or short-range correlation in momentum space: 

(qi,q2;Pi,P2) c ^ for |pi - qi| > l/\L\ and |pi - q 2 | > 1/|L|. (A. 9) 

In this case, using ( |A.9[ ) in ( |A.2| ), we see that ( |A.5| ) becomes, with obvious notation, 

5 °4 I ^^^ + ^-E^-E^^^^\^ 9Vi ^ (A.10) 

where V = L l L 2 L 3 . On the basis of the observation made at the footnote in §2.2.2, we 
assume that in flA.q* ) is insensitive to the region |qi| < 1/|L|. Then, in ( |A.6| ), the 
integration over qi yields an 0(T 3 ) quantity, with T the typical scale(s) of the system, 
which characterizes the microscopic reaction and is much larger than 1/|L|. Since we 
have assumed that A and B are of 0(1), we obtain 

S (n) yx'^T^v, 

Thus, S/S^ « 1 and the initial correlation C 4 may be ignored. There are, however, 
several cases, in which this observation does not apply: 

• Dimensionless quantity A in (|A.10|) is large. 

• T in (|A.10| ) is extremely large at q x ~ pi and/or qx ~ p 2 . 

• In massless theory (m = 0), some quantities are sensitive to the infrared region (cf. 
the footnote in §2.2.2). If S is ( |A.5| ) is such a quantity, more elaborate analysis is 
necessary to see if C 4 may be ignored or not. 
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Appendix B: Some properties of the propagator and 
the self-energy part 

Self-energy-part (E) inserted propagator is obtained by perturbatively computing the 
formula for the full propagator: 

G(x,y) = -i(T c U{x)i\yyS dAz& ^)), (B.l) 



where C int (= C — C ) is the interaction hat-Lagrangian. From this, we obtain Gn + 
G22 = G12 + S21, of which (|2.9| ) is a special case. Applying this to the single- E inserted 
propagator (Q — > A • £ • A) and using ( |2.9| ), we obtain Yll j=i ^ij( x >y) — 0- 
Taking complex-conjugate of ([B.1Q , we see that 



iGn(x,y) = [1^22(^,3;)]*, iGn(2i)(x,y) = [iGi2(2i)(y, 
Using this for G — > A and G — > A • E ■ A, we obtain 

-i(E(x, y)) n = H(S(j/, x)) 22 ]* , -i(E(x, y))i 2 ( 2 i) = &)) 12(21)] * • 

Fourier transformation onx-t/ yields 

(E(X;P)) n = -(E(X;P)); 2 , 

(E(X; P))i2 , (Epf ; P)) 2 i : pure imaginary 

Using this in (^33|) , we obtain E A (X; P) = P)]*. 

In the case of physical- iV scheme, £ c in §4 contributes to E, — V(x, y) G E(x,?/) (cf. 
( |4.57| )). This "additional" contribution does not invalidate the above properties. 



Appendix C: Representation of A c in terms of a vac- 
uum-theory kit 

Here obtain the expression for A c and the number density in terms of the quantities 
in vacuum theory. 
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Straightforward but a bit lengthy calculation using ( [A. 1| ) and its hermitian conjugate 
yields 

J d\x - y y p ^\^\y) + f (i/)0(x)) 
zA c °)(P) + iA«(X;P), 



*A C (X;P) 
zA(°)(P) 



(C.l) 



7T 



— [«J(po - P P ) + 6( Po + E p )] 
tip 

2n5(P 2 -m 2 ). 



iA^(X;P) = 2tt J d 3 q 



x [( a t( p+ ) a (p_)) ( j(p - (£ + + £_)/ 2 )e^+- s -) x °-q-x] 

+ (& t (-p + )6(-P-))^(Po + {E+ + E_)/2)e l ^ + - £; -) Xo+q - x ] 

+ -(a t (p+)fe t (-p-))5(p - (E+ - E_)/2)e i[{E++E - )x °-^ ] 
8 

+~<&(-P-)a(p + )}5(Po - (P + - £_)/ 2 ) e ^+*->*-**l 

'(C.2) 

where p-t = q ± q/2 and P± = P p± . 

Let us first show that the contribution of the third and fourth terms of the bracketed 
quantities in ( |C.2| ) is negligibly small. From the setup in §1, Ac (X;P) varies slowly 
in X, which means that (a^(p + )6^(— p_)) ~ (b(— p_)a(p+)) ~ for |q| > 1/|L| and 
E + + P_ > l/L°, where L M is the size of a spacetime cell (see §1). This can be realized 
only if m < 1/L°. However, as mentioned at the second footnote in §4.3, in the case 
of m < the bare propagator should be replaced with the tadpole self-energy-part 

resummed propagator. The tadpole induces the mass of O(VAT) with T the typical 
scale(s) of the system, so that E + + P_ >> 1/L°. 

Let us turn to the first and second terms in (|C2 ). Same reasoning as above leads 
to (a^(p + )a(p_)) ~ p+)&(— P-)) ~ for |q| > 1/|L|. We shall consider the "hard 
region" |P M | > 0(1/1/). From the first term, we pick out 



Pa 



-iq-X 



5{p - E p 



Ey - (p ■ q) 

8P3 



-6'(p -E p )+ ... 



-iq-X 



P 2 V 2 

8(p -E p ) + 8'(p -E p )-^- 



(P • V x ) 2 



8E* 



+ 



-iq-X 
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Then, to the gradient approximation, we have 5(po — (-E+ + P_)/2) ~ 5(po — E p ). 
Similarly, we have E + E_ ~ E'i. The second term in ( |(J.2| ) may be dealt with similarly. 
After all this, for |P M | > 0(1/L M ), we finally obtain 



tA c (X;P) ~ 2vr 0(p o ){l + 2ivt +) (X; p)} 



+0(-p o ){l + 2iVir ) (X;- P )} <J(P 



m 



(C.3) 



where 



< +) (*;p) 
a£°(*;-p) 



d 3 ?(a f (p + q/2)a(p - q/2)) exp 



/ < - q • X 



d 3 ?(& f (-P " q/2)6(-p + q/2)) exp 



Pr> 



/ < + q X 



E n 



are the number densities of the quasiparticle (+) and the anti-quasiparticle 
parison of O with (|2T22|) shows that 



(C.4) 



— ). Com- 



9 { b\X;p) 



N B + \X- P ) 



l-N { -\X;-p). 



For A C (X; P) with |P M | < one should use the fundamental formula ( p.2|) . 

Incidentally, in the case where p is diagonal in momentum space, (cr(p)a(q)) = 
5(p — q)n + (p) and (6t(p)6(q)) = 5(p — q)n_(p), we obtain iVg (X; p) = n±(p), where 
n.± are the number densities of the (anti)particle in vacuum theory. 



Appendix D Derivation of (|3.48|) and ( |3.50| ) 



Dl Derivation of ( ggg) 



Let us start with computing the piece of (j3.47T) , which diverges in the narrow-width limit, 
ImTiR — > — e(po)0 + . The relevant region (of integration) is where Re(G^ 1 (X; P)) = 
Re(Gf(X;P)) = P 2 - m 2 - ReY> R (X;P) ~ 0. We define "on the mass-shell" p = 
±uj±(X] ±p) (= ±u;±) through 



[P 2 -m 2 -PeS^(X;P))] 



= 0, (D.l) 

PO = ±OJ± 
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from which we obtain 



±Z ± 1 u±v± = p + 
dReZ R 



1 dReZ R 



2 dp 

2Zg 1 u± 



PO = ±U>± 



po = ±u>± dX^ 
Here v± = ±<9o;±/<9p are the velocity of the ± mode with momentum ±p and 

1 dReE R 



(D.2) 



Z7 1 = 1 T 



Po = ±w± 



(D.3) 



2^ ± dp 

are wave-function renormalization factors. To extract the "diverging" piece, we can 
make following approximations; 



F"£ o// ~ 0(po) K, pf sj} + 0(-po) (-"-, p)" 



-^[G R (X;P)-G A (X;P)]~^ 



#(rp ) ImE}j 



7^ [2cj t Z t ^po - ru; r )] + (IrriE 



where 



S o// (X;p = ±w±,p), = £ R (X;p = ±w±,p 



Using ( |D.4| ) in ( |3.47| ) , we obtain 



(SnX)) ~ - / 

W V ;/ 2J (2tt) 3 



(+) 



(-) 

off 



(D.4) 



(D.5) 



;d. 6 ) 



where V± = (1, ±p/a>±). In deriving (|D.6| ), we have assumed that | Jm£ jj| << 
under which ( p.4|) is a good approximation. Referring to ( |2.25|) , we can extract from 
( p.6| ) the contribution to the physical-number density: 

7 -?S (±) 

6N^\X;±u ± ,±p)- Z± ° ff 



2 l/mE^f 



D2 Derivation of ( g^OQ 

We note that r±V±/2 in ( p. 48 ) is the change in the "physical" -number density, during 



the time interval r±/2, due the the net production rate. In §1, we have introduced the 
spacetime cells, whose size is (// = 0,1,2,3). It is quite natural to take t± = L , 
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where Lq is the size of the time direction of the spacetime cell (including the spacetime 
point X^ in ( |3.48[ )). [Strictly speaking, in general, r + ^ r_. However, in the present 
crude argument, we ignore this difference.] It is interesting to note that r±V±/2 in 
( [3.481 ) is half of the net production probability, due to the reactions, during the time 
interval t± = L$. [In this respect, cf. Ref. 21).] 

Noticing that we are concerned about the quasiparticle mode with momentum p 
[anti-quasiparticle mode with momentum — p], we see that Q3.48 ) leads to the relations, 



8N±(X + L ,X + v±L ; ±w±, ±p) - 5N±(X , X; ±w±, ±p) 
-iE of f(X; ±w±,p) 



LnZ- 



2uo± 



from which we obtain 

d5N±(X;±uj ± ,±p) 

dX~ ' OX 

Using flDj) with Q, we get 



d5N± 
+ v ± • ~ 



1 ± + ' 



2u 



± 



-2Z± 1 cj ± 



po = ±w± 



p \ diV^pfjiu^ip) 



du± dNg\X;±po,±p) 
dX» dP„ 



± 



dX 



PO = ±0J± 



;d.7) 



(D.8) 



where use has been made of (|2.29 ) and ( |2.27 ). Substituting this into (|D.8j ) and adding 
( |2.29 ), we finally obtain 



dX 



+ v± • V x iV (±) 



' r{p) ^dcu ± dN { ±\X;±p ,±p) 



dX„ 



J P0=±U)± 



where = N ( ^\X; ±w±, ±p) + 5iV ± (X; ±w ± , ±p). 
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